NASA CONTRACTOR 
REPORT 


NASA CR-129008 


• H ^ 

tO O f 

a m 

as am 


m 

-O 

V 

ro 

CD 

rN 


AN EMPIRICAL ANALYSIS OF THE DISTRIBUTION 
OF OVERSHOOTS IN A STATIONARY 
GAUSSIAN STOCHASTIC PROCESS 


By Michael C. Carter and Michael W. Madison 
Department of Mathematical Sciences 
Appalachian State University 
Boone, North Carolina 28607 


July 1973 



in 

to n* u 
m < co • to 
to to & u 
as w 
^ h u 

< o u 

25 CO £33 3J 
<HOi 
C CL 

^ o u 

3C H ^ 

U CO H 00 
H £33 IT 
03 W < 

M t> 33 ^ 

Pm O CJ * 

E O > 

W CL. h H 

o to P 

25 » 

< as sc 

O < 0) 
HHL 
, — H co rtJ 

00 D CO 4J 
O CQ a? CO 

O M < 

o c 
NH ro 
»“* t/3 X *H 

1 H 03 

Q3 Q -fl) U 
U SC rtj 

* W O -4 

< a h <tj 
V) H H a, 

«c Cl 

* H < 
w O CO s *’ 


Prepared for 


NASA-GEORGE C. MARSHALL SPACE FLIGHT CENTER 
Marshall Space Flight Center, Alabama 35812 



CHAPTER 1 


INTRODUCTION 

The purpose of this thesis is to report the results of an 
empirical analysis of the frequency of overshoots above an 
arbitrary level in a stationary gaussian stochastic process. 

The problem is of interest to the Terrestrial Environment 
Branch, Aerospace Environment Division, Aero-Astrodynamics 
Laboratory, George C. Marshall Space Flight Center, Alabama, 
and the financial support for the project was under NASA 
contract no. NAS8-29286. The results obtained in this analysis 
are applicable in the prediction of extreme properties of 
processes such as wind speed, ambient temperature and sea 
state. The methodology may also be used in other fields, 
i.e,, electrical engineering and aerospace vehicle responses 
to forcing functions having known or assumed exponential 
autocorrelation functions . 

The remainder of this chapter presents a general dis- 
cussion of the scope of the work, and the organization of this 
analysis. 

1.1 Statement of the Problem 

The problem dealt with herein concerns the frequency dis- 
tribution of overshoots in a stationary gaussian stochastic 
process with an exponential autocorrelation function. Briefly, 
a stationary gaussian process may be described as a stochastic 
process which, at any point in time, has a gaussian distribu- 



tion. To be stationary the process must have a mean inde- 
pendent of time and an autocorrelation function dependent only 
on the distance between successive tine points. 

The problem has been of general theoretical interest for 
some time while little has been done to obtain numerical re- 
sults. Previous work in this general area is excellently 
summarized in two recent texts, Cramer' and Leadbetter (1967) 
and Kuznetsov (1965), and both contain extensive bibliogra- 
phies. The general density function for the number of 
crossings in a (0,t) time interval was given by Kuznetsov 
and Stratonevich (1956) . For a stationary gaussian process 
with R(t) = exp(-0x 2 ) Tikhonev (1956) approximated the 
probability of zero crossings in (0,t) by expanding the 
proof given by Kuznetsov and Stratonevich (1956) and neglect- 
ing terms in the series of order greater than 2. Other 
authors have various expressions for this density function 
and have investigated its asymptotic behavior. A general 
result states that as the level increases the number of over- 
shoots in (0,t) is Poisson distributed. A more extensive 
summary of previous work in this area is presented in Appen- 
dix II. To the author's knowledge this is the first investi- 
gation conducted by extensive simulation of such a process. 

1.2 Organization of the Analysis 

Chapter 2 is a discussion of the simulation model and 
assumptions concerning the model. Using the methods 
developed in Chapter 2, several simulations were run on an 
IBM-1130 computer. The results and analysis of the simula- 
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tions are presented in Chapter 3 along with the resultant 
distribution equations. 

The modus operandi for NASA to apply this solution to 
their specific problems concerning atmospheric variables is 
presented in the concluding chapter of this analysis. 
Appendix I contains a computer program to utilize the 
algorithm obtained in this investigation. 


CHAPTER 2 


Model and Simulation 

The first step in this development of a solution to the 
overshoot problem was to define a mathematical model of a 
stationary gaussian stochastic process with an exponential 
autocorrelation function. In developing the model, the 
following conditions were assumed: 

1) The sample process had a multivariate normal 
distribution. 

2) The process was strictly stationary, i.e., the 
autocorrelation function R(t^,t_.) = R(t) 

where t = It .-t. | . 

D i 1 

3) The expected value of a random variable X at 
time t was 0, i.e., E(X(t)) = 0 where E de- 
notes the expectation operator. 

4) The covariance matrix, denoted E, was symmetri- 
cal and positive definite. 

5) The autocorrelation function, denoted R(t), 

was exponential in nature, i.e., R(x) * EXP (— 3 | t | ) . 
The notation X(t) will denote a stochastic process satis- 
fying the above conditions. 

The process was considered over a time interval [0,99] 
and a sample realization consisted of 100 equally spaced 
sample points in the interval. This permitted some generality 
in the analysis whereas for a specific application the range 
of interest would be some [0,T] interval. In this case X(t) 
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would be sampled at t o ,t l' * * * ft 99 w ^ ere = ^To£p T ' 

with a corresponding modification of the autocorrelation 
parameter 3. The method of simulation was given by Odell 
(1971) and a summary of that technique is presented in the 
following discussion. 

Let X = (X(t Q ), Xft^, X(t 99 )) ' (" denotes matrix 

transposition) , then the covariance matrix is given by 
E = (a i j ) = E (X’X'*) so that for 0 < i, j < 99 

<? i j = E(X(t.)X(t j ) ) = R(t i ,t j ) = R(t) where x = l^-t^l. 

Thus it follows that £ is formed by evaluating R(x) for 
0 < x < 99 giving 


R(0) 

R(l) 

R(2) 

R (99) \ 

R (1) 

R (0) 

R(l) 

R (98) 

R(2) 

R (1) 

R(0) 

R (97) 

• 

• 

• 

• 

* 

• 

• 

* i 

R (99) 

R (98) 

• 

• / 

R (0) / 


By assumption X satisfies a multivariate normal distribu- 
tion with mean u * 0 and covariance matrix £, denoted 
X /^N(w_ r I). The following result, given by Odell (1971, 
pg. 37), provides the modus operand! of generating realiza- 
tion of X (t) . 

Theorem: If the 100 x 1 vector Y ~ N(u,E), and 

X is a fixed 100 x 1 vector, then 
V « AY + ]£ is distributed N(Aj^ + ^, A E A"). 
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We generated a vector Y ~H(0,I), I denoting the identity 
matrix, and obtained a factorization AA' of £, therefore, 
by the above theorem V = AY was distributed N(0,I). The 
resultant vector V constituted a realization of X(t). 

The generation of the vector (y fl ,y jf...»yg g ) > = Y^N(0, 1) 
was accomplished by generating a sequence of 100 independent 
standard normal variates. The Crout method was used to 
factor E into AA'. 

The technique of generating Y was given by Hamming 
(1962). He notes that an approximation to normally distri- 
buted random numbers can be produced from a sequence of 
uniformly distributed random numbers by the formula 

K k 
E *v " o 

y^ = k=l where x^ is a uniformly distributed 

Sk712 

random number in (0,1), and K is the number of values of 

x. used. According to the Central Limit Theorem, as K 

K 

tends to infinity the value of y^ approaches a standard 

normal distribution. To implement this procedure on a 

computer we fixed the value cf K at 12. The formula for 

12 

y. could then be expressed as y. * T, x v - 6. This con- 

1 1 k=l 

struction of y^ for 0 _< i _< 99 produced a sequence 

Yo'^1' ‘ * * '^99 stan( * ar< * normal variates with mean 0 and 

unit variance. This sequence is the vector Y ~N(0,I). 

We generated 250 realizations V^, i * 1,2,..., 250, 

for each of the autocorrelation functions simulated. This 

required 250 random vectors Y^, i » 1,2,..., 250 which in 
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turn, required a sequence of 250 * 100 = 25,000 standard 
normal variates, or 250 * (100 * 12) = 300,000 uniformly 
distributed random numbers. The algorithm used to generate 
uniformly distributed random variates was: 

r n = normalized (S n ) where S n = Xr n-1 

and the normalization is a reduction to (0,1). 

This algorithm is the well established power residual method 
of generating pseudo-random sequences. The period of the 
sequence generated in this fashion is a function of the in- 
teger capacity of the computer being used for the generation. 

In the case of the IBM-1130, the largest integer, and hence 
the period of the sequence, was 32,767 which falls far short 
of the necessary 300,000. 

Since the period of one number generator is too short 
to produce 250 realizations, we used a separate random 
number generator G ^ , j = 0, 1, ..., 99, for each of the 100 
elements of = (y A Q , y^, ..., y i ^ 99 > 1 < i < 250. Thus 

the generator Gy 0 _< j _< 99, produced the sequence 
Yl j' ^2 j' *•*' ^250 j ^dependent standard normal variates. 

In this fashion, each generator was required to produce 

250 * 12 * 3,000 uniformly distributed numbers, which is 
easily possible on the computer used in this analysis. 

To transform each vector Y^ into a realization 
of the process X(t) via the linear transformation * AY^ , 
it was necessary to factor the variance covariance matrix E. 
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As noted previously, £ is a symmetrical, positive definite 
matrix. A well known theorem in matrix theory states that 
such a matrix can he factored into the product of a lower 
triangular matrix and its transpose. This factorization, 

£ = AA' , where A is lower triangular, was accomplished 
using the Crout factorization technique as presented by Odell 
(1971, pg. 38). The method is summarized in the following 
discussion. 

The elements of A = (a^) will be computed in the 
following seouence: * * * , ^ 100 , l , ^ 22 , ^ 32 , * * * f ^100^ 

..♦»a 99 , 99 ,a 1 () 0 , 99 ,a 100 , ino . Note that A is lower tri- 
angular so a. . - 0 whenever j > i. Using this fact we have 


o . . 

ID 


3 


I a 
k=l 


ik 3 jk 


( 2 . 2 ) 


from which the following algorithms were derived. For 

2 

i = j = 1 we have = a^ so it follows that 


‘11 


(o 
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For i > j * 1 we have 


ij 


ments of the first column of 



(2.3) 


■ a ii a ii 80 ^e remaining ele- 
A are given by 


a il " °il /a ll- 

After j-1 columns of A have been generated we have 


(2.4) 


'jl 


k -i a j* a jk 


3; 1 2 2 

k-l** k + *11 


so for the remaining diagonal 


elements we have 


D 

For the remaining elements we have a . . = Z a. a . 

ID k=1 ik jk 

j-1 

£ a.,a.. + a..a.. so we can conclude 
k*l 1K ^k 13 33 

j“l 

= (0 ii “ 1 a nV a *iV )/a - • f ° r (2.6) 

iD 1 3 k=l lk DD 

i = j + 1, j + 2, . . . , 100. 

The autocorrelation function R(t) = exp(-B|T|) determines 
the degree of association between successive values of X(t). 
The process X(t) was simulated for a range of B values 
yielding processes where the correlation was above .98 
throughout the process, to processes where X(t) values 
could be considered independent after two time intervals. 

The minimum B value used was .002 which yielded R(99) = 
.9802, and the maximum B value was 5.0 which yielded 
R ( 2 ) = .00004539. The primary B values utilised were .002, 
.005, .0075, .01, .025, .05, .075, .1, .25, .5, .75, 1.0, 1.5, 

2.0, 3.0, and 5.0. We did, however simulate processes which 
were outside our primary range of interest, namely 7.5 and 

10.0. At each of these B values 250 realizations were 
generated. The selection of 250 as the number of realiza- 
tions for each B value was based on available computer 
storage capabilities but, from a statistical viewpoint, was 
deemed adequate for subsequent estimation and inference 
activities . 


v 
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CHAPTER 3 


Simulation Results and Analysis 

Once the "data sets" had been generated the basic 
problem of counting overshoots came into focus. Letting A 
denote some arbitrary level, we counted the number of over- 
shoots above values A = .5, .75, 1.0. 1.25, 1.5, 1.75, and 
2.0. Since each realization has mean 0 and unit variance, 
this was equivalent to counting the number of overshoots over 
.5 standard deviations above the mean, .75 standard devia- 
tions above the mean, etc. In future applications the over- 
shoots above a value of, say A - .75, would be equivalent to 
overshoots above a value of ,75a + u, where the process has 
mean u and variance o 2 . 

The value of 2.0 was selected as the upper limit of the 
major range of interest since, in the completely independent 
case, only 2.27% of the values would be above 2.0 and in the 
more correlated cases, the number of points, and hence the 
number of overshoots, would likely decrease. The value of .5 
was selected as the lower limit of the range of A values. 

In the completely independent case 30.85% of the values lie 
above .5, but the more memory the system has the longer the 
duration of each overshoot, and hence the fewer the number of 
overshoots. We did, however, count overshoots above higher 
levels for the purpose of determining the integrity of the 
estimation model outside the primary range of interest. 
Specifically, overshoots were counted for A levels of 2.25, 
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2.5, 2.75, 3.0, 3.5 and 4.0 for 0 values of .005, .02, .05, 
.1, .5, 1.5 and 3.0. 

To count the number of overshoots above level A, we 

counted the number of times V(t. ) < A while V(t.) > A 

l-l - i 

where 0 _< i _< 9$ and VCt^) ” 0* 

After the number of overshoots for a particular level A 
and autocorrelation parameter 8 was determined, the sample 
mean, X, and variance, S 2 , were computed in the traditional 
fashion. This provided the data to complete the table of 
means for A and 0, Table 1, and the table of variances 
for A and 0 , Table 2 . 





1 

* 


As noted in the introductory chapter, the Poisson dis- 
tribution is the limiting distribution as the crossing level 
becomes large and it seemed reasonable to first try the 
Poisson as a model for lower crossing levels. An estimation 
model for the multivariate Poisson (multivariate in the sense 
that the parameter X was assumed to be a function of A and 
3 ) was implemented and tried for various functions of A and 
3 . The results were discouraging. We first attributed the 
failure to our inability to find the proper function of 
A and $ , but later it was determined that the Poisson model 
was, in general, inadequate. 

The next and most fruitful step was the careful examina- 
tion of the means and variances for various levels of A and 
3. This led immediately to the following conclusions: 

1) There was a strong empirical relationship between 
the sample means and A and 3, and to a lesser 
extent, between the sample variances and A and 3. 

2) The binomial and negative binomial distributions, 
with parameters calculated from the sample means 
and variances, were more appropriate for the levels 
of A we investigated. 

For values of A < 1.5 and 3 ^ 1,0 the means exceeded the 
variances with the discrepancy increasing as A and 3 de- 
creased. As A and 3 increased above 1.5 and 1.0 respec- 
tively the values became approximately equal or the variances 
exceeded the mean. Once this trend was noticed, the reasons 
for observations 1 and 2 above became clear. If we assume 
one of three models, binomial, Poisson, or negative binomial, 
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is appropriate, then an accepted selection criterion is the 
relationship between the mean and variance. These observa- 
tions led us to seek those functional relationships that 
could best predict a process mean and variance. 

In the search for the relationship between A, 3 and 
the mean y, we first graphed the sample mean X as a 
function of A for each 3. This graph. Figure 1, strength- 
ened the conclusion that such a relationship existed but, due 
to our inability to find an appropriate approximating 
function of that relationship, this method of viewing the 
data was abandoned. However, we did note from this plot that 
the relationship behaved in what appeared to be an exponen- 
tial fashion. 

Suspecting the exponential characteristic, the next step 
was to graph ln(X) as a function of A for each 3 on semi 
log graph paper. This plot. Figure 2, was not a straight 
line as we had anticipated, but rather it seemed parabolic 
with the parabolas opening about the ln(X) axis. We selected 
the general parabolic model 

2 

In (X) = X Q (3) + X^gjA + X 2 (B)A (3.1) 

to try as an approximating relationship. The least squares 
technique summarized below was used to estimate X Q (3), 

X^(g), and X 2 (B) ^ or eac ^ Using these results, equation 

(3.1) was then rewritten to produce the estimate of the mean 
as 

EST(y) » exp (X q ($) + X 1 ($)A + X 2 ( 3 ) A 2 ) . (3.2) 
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For each 3 the estimate of y was a good approximation of 
X so we concluded that if the dependency of the A's on 6 
could be found, then (3.2) would provide a good estimate of 
the mean. 

The least scruares technicrue was given by Jorgenson 
(1961) and was used to estimate A = (Aq,A^, . . . , A k ) for the 
general model n^ = Agt^ + A 1 t 1 + ... + A^tj, for i=l,...,m 


(m being the number of observations). Let T = (T^) where 
= tj for j=l,...,k, and i=l,...,m. Then A = (T'T) - 1 T'n 
where n = (n, ,n_,...,n )' and A is a k x 1 vector of the 
estimates. 

Using the results of the least squares method, the first 
step toward determining the dependence of A^ on 3 was to 
plot In ( 3 ) vs. A^ on semi-log graph paper for each of 
Aq,A^, and A 2 . For all three coefficients three distinct 
trends were observed. For 3 ^ .01 the relationship was 
linear, for .01 < 3 <_ 1.5 the relationship appeared quadra- 
tic, and for 3 > 1.5 the relationship was again linear and 
essentially horizontal. Accordingly, the following models 
were fit using the method of least squares: 


3 < .01 In ( 3 ) = a Q + a^A i , i=l,2,3 

2 

.01 < 3 £ 1.5 ln( 3 ) = a + a A. + a A., i=l,2,3 

3 > 1.5 3 = a Q + a^, i=l, 2,3 (3.3) 

The results of the least squares estimates of the a^'s are 
presented in Table 3. 
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Using the results given in Table 3 we then solved (3.3) for 

Xq, X^, and yielding estimation equations for Xq, X^, 

X _ as 
2 


B _< .01 X i 

.01 < B <_ 1.5 X i 


B > 1.5 X. 


(In (B) - a Q )/a 1 i=l,2,3 

2 1/2 

(-a 1 + [a 1 -4a 2 (a Q -ln(B) ) ] >/ 2a 2 i:=1 ' 2 

(-a - [a 2 -4a (a -In (B) ) ] 1/2 ) /2a i=3 

112 0 2 


(e - a o )/a i 


i=l , 2,3 (3.4) 


We then estimated all coefficients and used them in the mecin 
prediction equation (3.2). The estimated means are given in 
Figure 1 along with the sample means. From this graph it is 
clear that, in almost all cases, the deviations are very 
slight and, as will be subsequently noted, the means estima- 
tion was deemed adequate. 

In the search for the functional relationship between 
A, B and the variance, a 2 , we graphed the sample variance, 
S 2 , as a function of ln(B) for each value of A. On careful 
examination of that graph, Figure 3, the ensuing observations 
were made; 

1) The sample variances were much more erratic 
than the sample means. 

2) For levels of A below 1.0, the graphs of the 
relationships of the variances and $ are, for 
all practical purposes, coincident (for this 
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reason only A = .75 vras graphed as a repre- 
sentation of all A <_ 1.0). 

3) The graphs are parabolic in appearance, opening 
about the In (6) axis. 

These observations led to the following model to estimate the 
relationship: 

In (3) = Aq (A) + A 1 (A)S 2 + A 2 (A)(S 2 ) 2 , where 
S 2 is the sample variance. (3.7) 

We ran least squares fits for each A and found that 
the model was acceptable provided the dependency of A Q , A^ , 
and A 2 on A could be determined. Toward that end, we 
graphed each coefficient of (3.7) as a function of A. The 
plots of Aq and A^ appeared linear. The graph of a 2 
at first appeared to be quadratic, but was later found to be 
better approximated by a cubic equation. Therefore the 
following models for A Q » A^ and A 2 were fitted using least 
squares techniques: 

x o * \ + a i A 

A, * a A + a. A 
10 1 

i 2 3 

' X 2 = a 0 + a i A + ®2 A + a 3 A * (3.8) 

The results of the least squares fits provided the following 

estimation eauations for A n , A, , and A : 

u ± 2 

A„ - -7.013 + . 3871A 
0 

\ 1 • .2192 - J.759A 

A 2 - -6.1371 + 14.813A - 11.633A 2 + 3.05A 3 . (3.9) 
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1 

.... . . 4 

Tests of (3.9) yielded good approximations to the coefficients 
for each level of A. 

The original variance model (3.7) was then re-written as 

0 = (X Q - ln(3) ) + XjS 2 + X 2 (S 2 ) 2 (3.10) 

and solving (3.10) for the estimated variance o 2 , we have 

EST. o 2 = (-X + [X? - 4X (X . - ln(8))l 1/2 )/2X . (3.11) 

112 0 2 

Using the coefficient model to estimate X , X. , X , we 

0 1 2 

then estimated the variances from (3.11). As noted previously, 
the sample variances are more erratic than the sample means 
and the estimated variances were not, in general, as accurate 
as the estimates of the means. The estimated variances are 
shown in Figure 3 along with the sample variances. 
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To check the accuracy of our results, we ran estimations 
of overshoots for the values of A and 0 that were used in 
the estimation process. The estimated mean was calculated 
from equation (3.2) while the estimate of the variance was 
given by (3.11). To determine which distribution was appro- 
priate, we formed the ratio 

e stimate of mean . . 

r “ estimate of variance* (3.12) 

If r <_ .95 then the variance clearly exceeded the mean so 
we used the negative binomial distribution 

P(X=i) * rfcf ri pkqi; OiPi 1 ' p +< 3* 1 » k>0 ' i-0,1,2...- (3.13) 

If .95 < r < 1.05 then the mean and variance were approxi- 
mately equal so we used the Poisson distribution 

*V X 

P(X-i) - “ — X>0, i«0,l,2,.... (3.14) 

If r ^ 1.05 then the mean clearly exceeded the variance so 
we used the binomial distribution 

P(X»i) ■ (^ip^q** S P+q-1* i«0,l,»**,n. (3.15) 

The test used to check the goodness of fit for the pre- 
dicted models was the Kolmogorov goodness of fit test. 

Briefly, the test compares the theoretical and sample distri- 
bution functions and one concludes there is no significant 
difference between these distributions if the maximum absolute 
difference between them is less than a predetermined quantity 
based on the significance level and sample sise. 
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The test is designed to compare continuous distribution 
functions and, as such, is not direct 1 ' applicable to dis- 
crete es. When it is applied the si gnificance level used 
is conservative to an indeterminate dc,,< -se . For our purposes 
this is quite acceptable. In g eiier.s. a conservative test 
conducted at a level of signi r ianr is, in reality, being 
conducted at some a" < a lew ; of significance. Thus in 
Table 4 the a levels are given as a < .05 or a £ .01. 

This means that a conclusion that we have a good fit using 
the o * .05 significance level really says the two distri- 
bution functions are in agreement at some a" value smaller 
than .05. For the cases in Table 4 where the predicted 
models fitted poorly we can only state that the model was 
rejected at some a level less than .01. The justification 
for using the Kolmogorov test for these data is given in 
Noether (1967, pp. 17-18). 

The sampling distribution for the Kolmogorov test is 

well known and for sample sizes above 35 the maximum absolute 

difference between the theoretical and observed distribution 

functions must not exceed the value d //n where 4 - 1.36 

o .05 

and d * 1.63. These results are available in Siegel 
(1956, pg. 251) . Consequently using n - 250 the critical 
values are 21.50 and 25.77 respectively. 

A brief glance at Table 4 shows excellent results through 
most of the 0 and A values with no rejections in the 
.1 < B < 1.0 range which will be the primary 0 values used 
in wind speed calculations. It is unlikely that many appli- 
cations will require those 0 values giving poor results, 



25 



namely .002 and 2.0. The .002 data set is the "end point" in 
our predictive process and 2.0 is the data set just above the 
"transitional" value where the means behavior changed 
drastically (previously discussed in this chapter). 

Appendix III presents a spectrum of data sets in the 
computer format used to evaluate the goodness of fit. All 
pertinent information, i.e., A and $ values, observed and 
predicted means and variances, model utilized, cumulative 
distribution functions and predicted probabilities, are 
presented. 
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TABLE 4 


Condensed Summary of Fitted Models 
Using Predictive Equations 

♦Good Fit at a £ .01 y + Bad Fit due to mean overesti- 

mation 

♦♦Good Fit at a £ .05 

U+ Bad Fit due to mean under- 
estimation 

a 2 + Analogous to definitions on 
a 2 4- y above 



.5 

.75 

1.0 

1.25 

1.5 

1.75 

2.0 

2.5 

3.0 

.002 

** 

** 

** 

** 

y + 

y + 




.005 

** 

** 

** 

* 

** 

** 

* 

y + 

** 

.0075 

a 2 + 

if it 

** 

** 

** 

** 

* ★ 



.01 

a 2 + 

if 

if if 

** 

** 

** 

** 

** 

** 

.025 

** 

if 

* * 

* 

* 

** 

* 



.05 

** 

if if 

** 

** 

** 

** 

** 

y + 

★ * 

.075 

** 

** 

** 

** 

** 

* 

* 



.1 

** 

if if 

if* 

** 

** 

** 

** 

** 

** 

.25 

* 

** 

** 

** 

** 

** 

** 



O 

u> 

• 

** 

** 

** 

* 

* 

** 

** 

** 

** 

.75 

** 

** 

** 

** 

** 

** 

** 



1.0 

** 

if if 

** 

** 

* 

** 

** 



1.5 

** 

** 

** 

** 

** 

** 

** 

** 

** 

2.0 

y+ 

y + 

* 

* 

* 

y + 

U + 



3.0 

#* 

** 

** 

** 

** 

** 

** 

** 

** 

5.0 

y + 

V + 


** 

** 

** 

** 
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The fact that the binomial and negative binomial models 
fitted the data was, upon reflection, not surprising for the 
following reasons: 

1) A well known statistical fact given by Johnson 
and Kotz (1969, pg. 43) states that, if one of 
the models is applicable, the criterion for 
selection depends on the relationship between 

the mean and variance, i.e., if y > a 2 select the 
binomial, if y is approximately equal o 2 
select the Poisson model, if y < a 2 select the 
negative binomial. 

2) As pointed out by Johnson and Kotz (1969, pg. 135) 
for 6 values that are quite low the negative 
binomial is the appropriate model since, with a 
small g value, the successive time points and, 
therefore, successive overshoots are dependent. 

In applications where the Poisson model seems 
appropriate but successive events are not indepen- 
dent the negative binomial model is an excellent 
alternative. 

3) For larger g values the binomial model is re- 
quired since successive time points and over- 
shoots are, for all practical purposes, independent. 

The criterion for a successful model was adequate fits 
on the majority of the data sets. As pointed out in Table 4, 
the experimental results for means and variances were, in 
general, approximated adeauatelv by the prediction model. 
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Therefore we concluded that the model presented in this 
analysis is a good predictor of overshoots in a stationary 
gaussian stochastic process with an exponential autocorrelation 
function. 


29 



CHAPTER 4 


Applications 

The purpose of this chapter is to explain how the 
analyses discussed previously could be applied to problems 
of a general nature where the assumption of a stationary 
Gaussian process with an exponential autocorrelation function 
is plausible. 

As noted in chapter 3 the A levels of .5, .75, 1.0, 
1.25, 1.5, 1.75 and 2.0 were used in obtaining the basic 
prediction models. Figure 1 gives the predicted means for 
levels above 2.0 for representation values of £. It is 
apparent that the predictive equations are adequate for A 
levels above 2.0. 

While the study involved counting overshoots above 
specified A levels it is valid to assume the model is 
applicable to predicting the number of "undershoots” below 
negative A levels - for no reason other than the symmetry 
of the normal distribution. 

The estimating equations for the mear and variance 
were derived based on realizations of 100 points. Some 
applications, based on time periods yielding appreciably 
different numbers of points in a realization, require a modi- 
fication in these estimating equations. In the computer pro- 
gram this is cone directly using the standard formulas. 
Letting M be the number of time points desired we have 
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New EST(y) - E;T(u)*M/100 


and 

2 

New EST(0 2 ) = EST (a 2 ) * (M/100) 

with the model being selected based on the new estimated 
mean and variance. It should be noted that these new 
estimates are exactly correct only for the independent case 
(high A and/or 8 values'. When the correlations go to 
zero rapidly (8 values >_ 1 and/or moderate A value*?' 
the new estimates will be higher by an unimportant, and 
likely indeterminable, amount. For small 8 values the new 
estimates could be inflated if A is also low. As our 
primary range of interest has been in 8 values _> .1 and 
A values distant from the mean the modification in the 
mean and variance will be satisfactory. 

The computer program given in Appendix I has been de- 
veloped to support applications of this study. We will, in 
the ensuing discussions, relate applications that can be 
performed using this program. 

To utilize this program the following data must be pro- 
vided : 

1) the average, y, of the process, 

2) the standard deviation, a, of the process, 

3) the coefficient, 8* of the autocorrelation function 
R(t) - exp (-8 1 t | ) , 

4) the number of time points and 

5) the crossing level L and maximum frequency N. 
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Assume that the probability distribution of the number of 
overshoots above some level L > y is desired. The pro- 
gram will use L, y, and a to calculate the A level, 
i.e., A = (L - yj/o, which will be used in the prediction. 

It is important to note that the model expects A to be 
positive and should L be less than y the A value would 
be | L - y | /o . 

The $ value used will, in most cases, correspond directly 
to the 6 values used in the analysis since empirically g 
values are calculated or estimated using serial correlations 
of lag 1, lag 2, etc. which are independent of the interval 
between successive time points. Shi la g value be calcu- 
lated using the actual time intervals it will be necessary 
for the user to modify the g value prior to utilizing the 
program. Recall that the g used in the program assumed "time" 
units of length 1. If a g value has been calculated using 
intervals of, say, .5, i.e., r = .5, 1.0, 1.5, etc. the auto- 
correlation function will be R(x) = exp(-g|x|), x = .5, 1.0,... 
and this corresponds directly to R(x') = exp (— . 5g | x “ | ) , 
x >- l»2, ...* In this case the value ,5g would be the 
value the user supplies to the program. 

In general we can summarize this procedure as follows: 
Assume the autocorrelation parameter 0“ has been calculated 
using equally spaced intervals » h, 2h, 3h, . .., giving 
Rtx*) * exp(-$' , |x < ’| ) . This corresponds directly to 
R(x) * expt-g'hlx | ) , x « 1, 2, ... which means 0 (for 
program input) = g -h . 
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The program output will consist of: 

1) L, y , 0 , 3 and calculated A value, 

2) Predicted mean and variance for the number of 
crossings and the model selected based on these 
values and 

3) Predicted probabilities for 0, 1, 2, . .., 
n (or more), n < 40, overshoots. 

As an example consider the situation below. For the 
month of January at 12 km the scalar wind speed at Cape 
Kennedy has the following properties: 

1) K (t ”* ) = expC-B^T") , x" = 12, 24, 36, ... 
with 6' = .0247, 

2) a = 8 m/sec and 

3) y = 24 m/s. 

We desire to predict the probabilities of 0, 1, 2, 3, 4 and 
5 (or more) overshoots above the level L = 39 m/s. 

The 3 value for the program is not .0247 but rather 
is 12 * .0247 = .2964. This makes R(t') = exp (-. 0247 t ') , 
t“ ■ 12, 24, 36, ... equal to R(t) = exp(-.2964t) , 
t ■ 1, 2, 3, ...* The program input is L * 39, y = 24, 
a » 8, B ■ .2964 N ■ 5 and M * 62. The program calculates 
the standardized crossing level A as 1.875 and utilizes 
these A and 8 values to calculate the predicted probabili- 
ties. Table 5 gives the resultant computer output. 

The formats for program input parameters are given in 
Appendix I. 
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TABLE 5 


LEVEL 39.0000 MEAN 29.0000 STD. DEVIATION 
ADJUSTED LEVEL 1.87SO DlST. MEAN 1*1259 DlST. V AK 
AUTOCORRELATION PARAMETER U*2969 

MEAN EXCEEDES *ARlANCt> BINOMIAL MODEL SELECTED 


NUMBEK OF CROSSINGS 
0 

1 

2 
3 
9 
S 


PREDICTED PROBABILITY 
D.3001 
0.3898 
0.221S 
0*0719 
0*019* 

0*0021 


b.OOOU 
U • 9686 
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This algorithm was, of course, developed under the 
assumption that the process was stationary normal with an 
exponential autocorrelation function. This process is the 
most widely applied and, despite the intractability of the 
mathematics, in the only such continuous stochastic process 
that can be analyzed to any appreciable extent. The question 
of the degree of applicability to suspected or confirmed non- 
normal processes will certainly arise and while there appears 
to be no answers in the literature there are soma statistical 
results of a general nature that have some bearing. 

Most skewed distributions, e.g., lognormal or gamma are, 
for certain ranges of parameter values, almost normal and 
certainly the algorithm is useful in these cases. As an 
example of this term "almost" normal consider a gamma distri- 
bution with parameters a = 2f and 6 * 1/2 where 

0 a x ct ‘ 1 e“ Bx 

f(x) * TT. • As f increases this distribution 

T (a) 

rapidly approaches normality. With these values of a and 6 
this is the x 2 pdf and it is common statistical practice 
to use a normal approximation to determine critical values 
when the degrees of freedom (denoted by f) is large. 

The equations to estimate the process mean and variance 
do not utilize the properties of any distribution. The 
assumptions of stationarity and exponential-like autocorrela- 
tion function are certainly required whether the distribution 
is or is not highly skewed. The Central Limit Theorem (CLT) 
permits one to determine the average number of points above 
a certain level in this process just as if it were normal 
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f 


(the CLT is certainly applicable with sample sizes of 250) and 
an extension to the average number of overshoots is reasonable. 
In skewed distributions the mean and variance are not inde- 
pendent - this possibly makes our independent equations for 
estimating the process mean and variance a bit tenuous. The 
relationship would certainly be difficult to determine but 
the proper procedure would likely be to incorporate the esti- 
mated mean (which should be adequate if the assumptions are 
satisfied) into the estimating equation for the variance. The 
three discrete distributions, i.e., binomial, negative binomial 
and Poisson, are applicable in any case. 

One additional result is the fact that overshoots fre- 
quencies approach the Poisson distribution as the crossing 
level increases regardless of the process distribution. This 
"cutoff" value was approximately two standard deviations above 
the mean in our study and would likely be close to the "cutoff" 
value for most distributions. 
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APPENDIX I 

PREDICTION PROGRAM LISTING AND 
INPUT FORMATS 
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SOURCE - STATEMENT PA4E OU| 

< 

c 

c 

C PROGRAM 10* 

C MMU1* 

c 

C AUThOR* 

C MlAE MAlISOn* 

c 

C date AftiTTEN. 

C 0**/|U/73. 

c 

C PURPOSE* 

c calculate frequency oistmiluTIOn for overshoots p 

C A STATIONARY 6ASUSSIAN STOCHASTIC PROCESS «ITh 

c eaponential autoconhelaTIOn function, 

c 

C SUbHOUTlNCS REOUlHED. 

c oln«*: calculate the loo of the oahha oistrioution. 

C 

c 



c 

c 

c 

INTtfcE* ror.prntr .ratio 
real hean, level, JP l 

C SET UP 10 ASSI6NHENTS 

RDRa | 

PRNTR • 3 

C 

C REAO PARAMETERS AND CONVERT To STANDARD CROSS I NS LEVEL 

C 

1«0 REAO t ROM , 1 000 ,E NO* 7 00 1 HEAN.STdEV , RET A, LEVEL, NUN, RAT 10 
NUN • NUN ♦ I 

• ■ -beta 

APRT «ILEVEL • MEAN) / STDEV 
A • ABS I APRT I 
ALOCB • AL06 I BETA I 
C 

C COMPUTE EST OF DISTRIBUTION MEAN 

C 

IF (BETA *6T. 1.6) 60 TO 260 
IF tBCTA *6T * *011 60 TO 200 
FI a UL06B ♦ 6.0717) / I. IIOS 
F2 a IAL06B ♦ 4.46644) / 1-3*23)16) 

F3 • UL06B ♦ 7.41706) / (-i0.|T2BI 
60 TO 300 

200 01 a l*67464**2 - 4.0* *21 B6I • I -4*623 - AL06M ) 

02 • 10*31 IB4**2 - 4 *0*10 *23467* <- l *0602 - AL06B) 

03 * 3*302062**2 - 4.0*14. 3464*1-3. 4H2 - AL06B) 

FI « 1-1*67604 * SORT 10 1 ) ) / *43702 

F2 • 1-10*31164 ♦ SORT (02)) / 20.47314 
F 3 a 1-3*302062 - S0)TI03)) / 32*73)2 
60 TO 300 


./ 
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r\ n n o r» n 


SOURCE 


STATtMtNT 


f AGE 0U2 


2<i0 P| ■ 3.305 ♦ .U19.B 

F2 m -.lb? - • 0 1 Ob*b 

F 3 ■ - • M 7 ♦ • 02» b 
3 q0 temp « FI ♦ F 2 • A ♦ F 3 • A « • 2 
tMtAN ■ tXP (TEMP) 

EMEaN«LMLAN.KATIO/|uQ. 

COMpUTt tST OF DISTRIBUTION VAkIANCL 

IF (A .fit. 1.0) 00 TO 3SU 

G| = -6.32799 

02 - -.3902259 
G3 * >16156 
00 TO 9Q0 

360 01 - -7.1)12996 ♦ ,387066»a 

12 » .219193 - . 1 758552* A 

03 " -6.137068 ♦ 19.81297‘A - ) 1 . 632 VS* A ♦ *2 ♦ 3 . 0996 1 • A • »3 

9(lO DISC « 02**2 -9.0»03«(Ol - aLOuu) 

EVA* «(-G2 ♦ SU h 1 (DISC) ) / (2*03) 

EVAr«EVAW.(KAT10/10U.)»*2 

DETERN I Nt APPROPRIATE MODEL AND REASON AND PR T Fit AD I NOS 

RLTn » EMEAN / tV AR 
IF (A .GT. 2.0) 00 TO SOU 

IF (.95 .LT. RLTN .AND. RLTN .LI. 1.05) 00 TO bid 
IF (EVAR .OT. EMEAN) 00 TO 520 
C MEAN > VAR -> BINOMIAL 

MODEL • 1 

PP • l - EVAR / tMEAN 

XTEMP ■ EMEAN««2/ (EMEAN-EVAR ) 

XTEMP ■ XTEMP ♦ .5 
ITEMP ■ XTEMP 

IF (ITEMP .LT. (NUM-ll) 00 TO 990 
PX « ITEMP 
00 TO 995 
99O PX a NUM - 1 
995 CONTINUE 
00 TO 550 

C A > 2.0 •> POISSON 

SoO MODEL * 2 
00 TO 5SU 

C MEAN • VAR -> POISSON 

5 | 0 MODEL « 3 
00 TO 550 

C VAR > MEAN -> NEfa B 1 NOM I At 

520 MODEL ■ 9 

PK ■ EMEAN*»2 / (EVAR - EMEAN) 

PP a PX ( (PK ♦ tMEAN) 

550 WRITE (PRNTr»|0|0)LEVEl»MEAN*SToEV iAPRT .EmEAN.EVAR.BETA 
00 TO (551 .552.S53.5S9) .MODtL 

551 WRITE ( PRNTR .|0||) 

GO TC 560 

552 WRITE (PKNTR.I012) 
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i 


0 


n n r- non non nno 


SOUrCL 


statement 


PAGE OUj 


GO TO S6U 

*R1TE l PHNT R . 1 0 1 3 ) 

GO TO S60 

55<t WRITE ( PKNTK • I 0 1 *4 > 

S*0 WRITL (PnNIK.IQIS) 

C 

C CALCULATE AnD PRINT ALL BuT LAST CELL 

C 

PSUM * 0 

NM J ■ NUM - I 

IF (NUM .EG. 1) GO TO 60U 

00 £90 l ■ l * NM 1 

J » 1 - I 

JPi ■ J ♦ 1 

GO TO ($71 >S72»S73.S7H) .MODEL 
C BINOMIAL MODEL 

67 1 PXPi ■ P* ♦ l 
OMPP • I - PP 

ALPR • DLGGMIPXPl ).J.A10G(PP).(PX-J)*AL0G(0MPP)-DLGGM( JPI I 
1-DLGGM(PX-J*1.» 

GO TO beo 

C POISSON MODEL 

672 CONTINUE 

573 ALPR • J * ALOG(EMEAN) - EMtAN . DLGGh (JPI) 

GO TO SBO 

C NEGATIVE BINOMIAL 

S7H PKPj ■ PK ♦ J 
OMPP • 1 - PP 

ALPR » DLGGM(PKPJ)^PK*AL0G(PP)+J*AL06(0hPP)-0LGGM(PK)-DLGGM( JPI) 

CALC PKOBABILITy AND PRT CELL 

S«0 PROB ■ EXP (ALPR) 

PSUM • PSUM ♦ PR0B 
S 7 O WRITE (PRNTR.1020) JiPROB 

COMPUTE LAST CELL 

*00 PROB • l - PSUM 

WRITE (PRNTR. 1020) NM 1 .PROB 
GO TO 100 

END OF JOB 

OO STOP 


formats 


)000 FORMAT (rF10.3i2IS) 

10)0 FORMAT »*!*, 'LEVEL*. Fll.R*SX.*MEAN*,FU.**,|QX.*STD. DEVIATION** 
|F I !**♦/* ADJUSTED LEVEL* *P 1 1 .*• OX . *DIST. MEAN* *F 1 1 .*) ,§X , 

2'DIST* VAR**Fll.H/* AUTOCORRELATION PARAMETER «.Fll.«*/) 
toil FORMAT ( MEAN EXCEEDES VARIANCE* BINOMIAL MODEL SELECTED*) 

1 0 1 2 FORMAT (« LEVEL ABOVE 2.0* POISSON MOOEL SELECTED* I 


f 


I 


■w— 
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SOURCE 


STATEMENT 


PA6E 00** 


I u 1 3 FORMA I (• MtAN APPhOX a VARIANCE* RuISSON hoOEL StLECTl. D‘ ) 

1 0 1 H FORMAT (• VARIANCE l.ACEEDES MCAn, NE&AI1VE BINOMIAL MONEL SELECT** 
1 • ED * ) 

1 0 1 b FORMAT (/^X, ‘NUMBER OF CRuSSlNui PREDICTED PRO'S AB I L 1 T T • ) 

1020 FOKmAI (• * . 10X , 1 6 .21 X ,F7.R) 


END 




i 


i 




source 


STATEMENT 


page ooi 


FUNCTION OUoGMCUA) 

0T»0* 

DTtRM* 1 , 

IFtOAl 1 • 1 |2 

1 ULGGH»U. 

HE T Orth 

2 lFloT-18. 1 J.3*** 

3 dterm»oterh*ov 
dy*dy*i. 

GO TO 2 

** 0L6gM» 1DY-»*»>*AL0G (UTI“DY»l,/(i i ;.«Uf)-l,/(J60,»L)Y»»3)*i,/(i260.» 
|DY««SI~1 •/( 1680.*uT**7 14..91 tt»3bbJ32uHA73-ALO5(0TEKM) 

RETURN 

END 
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APPENDIX II 


COMMENTS ON THE THEORETICAL APPROACH 


COMMENTS ON THE THEORETICAL APPROACH 


Two interesting problems in the theory of stochastic 
processes are first to find the probability density of the 
duration of a crossing of a given level by a random process 
X(t) and second to find the probability density of the 
number of crossings of a level by the process. The problem 
of obtaining the average number of crossings of a level has 
received much attention in the literature. In fact, if 
X(t) is a stationary gaussian process, the complete solu- 
tion has been given by Ito (1964) and Ylvisaker (1965). For 
non-stationary gaussian processes, Leadbetter and Cryer (1965) 
have given a similar result. And finally, Leadbetter (1966) 
has considered the average number of crossings for a wide 
class of non-gaussian processes. However, solutions in 
closed form for the original two problems have not been ob- 
tained even in the more desirable case when X(t) is gaussian. 
Several approximations to these probabilities have been ob- 
tained and we shall give some with references. 

Let X(t) be a random process with correlation function 
R(t). Following Rice (1945), the probability density function 
of the interval between the ith and the (i + m + l)th cross- 
ing of a level A by X(t) is denoted by P (t); and the 

m 

probability of exactly n crossings of the level A in the 
interval (t, t + t) is denoted by p(n,x). For a basic re- 
lationship between p m ( T ) and p(n,T) see Appendix I of 
McFadden (1958). Let f r = X(t Q ) = A, ^ = X(t i ), gj^ = 
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for i = 0, 1, . .., n and d (t. , . .., t ) = 

K K 


(1/N (A) ) dg 


0 V J -«o 


^ g 0 g i*" g kl W k+i (A '“*' A 'g () ^-- c ? k )^g i *..dg k 


where N(A) = g W (A;g)dg is the expected number of cross- 
JQ 1 ' 


ings of the level A by X(t) in the interval (t, t + t) 

and W (f ,...,f ;g^,...,g, ) is the joint probability 
k+1 0 k 0 k 

density of fg,...,f k and 9o'***' g k* T ^en according to 
Kuznetsov and Stratonovich (1956) 


00 k f T f T 

p (n,x) = 1/ni E (-1) /k! ... d (t, , . . . t )dt . . .dt J _. . (1) 


0 JO 


n+k 1 


n+k 1 n+k 


And by Kuznetsov, Stratonovich and Tikhonov (1954) the 
probability density for the duration x is 


p(x) = p(0, T ) . 

dr 

It is apparent that the desired probabilities (1) and 
(2) are very complicated and conseguently only approximations 
have been given. To illustrate this point let us consider 
the case where X(t) is gaussian with correlation function 
R(x) . It is known that the normal property is retained for 
any linear transformation of a normal random function. Con- 
sequently, the joint probability density for the values of 
the random function and its derivatives will also be normal. 


Thus w k.i (f o f k !9 o V 


d/un) k+1 4 1/J 


2k+l 


) exp (-1/2 E L. . f , f . ] 
i,j*0 13 1 3 



I 


f = g , , 

k+i l-l 

i=l, . . 

• f k+l , 

A = | (r. .) | 
ID 

• • • 

00 

R 

0k 

• 

R' 

00 

• 

R" 

0k 

• 

o 

• 

n 

V. 


• 

• 

• • • R. . 

kk 

' " R nn 
00 

i 

! 

\ 

-R ai 

0k 

• 

-R'' 

00 

• 

-R' ' 
0k 

• 

\ 

\" R k0 * ' • 

-*L 

• 

-R'' 

k0 

• 

••• -*££ 


2k+2,2k+2 


and R. . = R(t . -t . ) . 
ID 1 3 


Thus we have W^ +1 (A, . . . , A;g Q , . . . ,g^) 


= (l/(2II) k+1 A 1/2 ) exp [-1/2 ( E A 2 L + E L v+1 + .g,gJ] 

i , j=0 i: > i , j=0 k+l+i f k+l+3 i 3 


2 k k 

= (1/(211) k+ 1 A 1/2 ) exp[-A /2 E L.J exp [-1/2 E L. . 

i,j=0 i,j=0 K+1+1 ' k+ 

Now let us denote by p (g Q , . . . ,g k | f Q , . . . , fy) the con- 


£ L g q ] 

i# j.o k+i+i,k+i+jVj J ‘ 


ditlonal probability density of (g^ , . . . ,g^) given (f , 
and let p(f , ...,f ) be the probability density of 

w Jv 

C f 0 # then 

p(f_,...,f t ) = (l/(2n) k+1/ V /2 ) exp (-1/2 E M . , f f , ) 

0 k i,j-0 L 3 

where (M ) is the inverse of (R^J and 


R. . • • • 

00 Ok 
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The last summation of (3) involves only the last k+1 rows 
and columns of (L ) . Denote the inverse of this matrix 

ij 

by (m^); that is, 


k+ 1 , k+1 


L k+l,2k+l \ 


-1 


(m. . ) = 
ID 


\ L 


\ 2k+l ,k+l 


J 2k+l,2k+l ' 


It is clear that (m ) is the covariance matrix of 

ij 

(g 0 ,...,g k ) given that f Q = f 1 - . . . - f k - 0 and by 

Jacobi’s Theorem the (i,j)th element of this matrix is the 
bordered determinant 


■« 


00 

• • • R 

Ok 

r' 

0j 

• 

• 

• 

• 

• 

• 

• 

• 

• 

‘k0 

R 

kk 

R' 

kj 

► 

-Rf 

-R'' 

10 

ik 

ID 


* D. 


The determinant of (m^) is given by I ( m i j > I = A / D * So 
now we have p (gQ» • • • #9^ 1 0, . . . ,0) 
y + l 1/2 ^ 

_ a/m) s ) «xpi-V2 i E_ n L lt+1+iiK+1+j g i g j ] 

, k+1/2 1/2 

(1/(211) n ) 

- (l/(2n) k+1/2 |(m lj )| 1/2 ) «Pl-W l( i 0 W l w + jVj l 

« z (g , m) . 
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V' 


i 


Thus W (A,...,A;g , . . . ,g ) 

k+1 0 k 


k +1 1/2 2 k ^ 

= (l/(2n) A / ) exp [-A /2 I L. .] exp[-l/2 I L . 

i,j=0 ^ i, j=o k+i+i,k+i+] i j 

k+1/2 1/2 2 k 

= (l/(2n) D ) exp [ -A / 2 E L.J . Z (g, m) . 

i,1=0 13 

Therefore, (t^ , . . . , t^) 

= (1/M (A) (2n) k+1/2 D 1/2 ) exp[-A 2 /2 E L.J.fdgf ...f |g 0 ...g| 

i f j=0 Jo U J-00 J — 00 K 

.Z (g, m)dg^. . .dg^. 


Let 


‘ij 


m ) 1/2 

ii n jj ) 

i + r. 

m ) 1/2 

ii j 

i = 0 


1/2 


and h^ « ^/(ro^)" then d^ (t^, . . . ,t^) 


((m 00 ,,,m kk )1/2/N(A) (2n)k+1/2r)1/2) exp[-A 2 /2 i E ^L i;j ] 


• | « • • | hg* • 'hjj^ n)dhg. . .dh^ 




where Z(h, n) is the ordinary normal probability density 
function in k+1 variables h n ,...,h k with covariance matrix 

<V- 

In an attempt to find p(t) when X(t) is gaussian and 
R(t) ■ exp(-ax 1 )# Tikhonov (1956) has approximated p(0,t) 
in (1) by neglecting all terms of the series greater than 2. 

He claims that his results give satisfactory agreement with 
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the experimental results of Rice (1953) . It is clear that the 

smaller t is, the better the approximation. However, if t 

is very large or if we v/ish to find p(n,t) for large n we 

must find some other means of approximation. 

Longuet-Higgins (1962) has obtained an infinite series 

for p(n,r) and P (t) similar to (1) where each term is 

m 

an integral of the joint probability W(+, -, -, ..., -)dt^...dt n 

that X (t) has an up-crossing in the infinitesimal interval 

(t^ , t^ + dt^) an< * a dow n-crossing in the remaining (n-1) 

intervals (t^, t^ + dt^) (i = 2, 3, . . . , n) . He -Iso gives 

a general relation between P m (T), p(n,r) and W(S) where 

S is a series of plus and minus signs (plus if X(t) has 

an up-crossing and minus if X(t) has a down-crossing). 

Using the infinite series he obtains the asymptotic behavior 

of P_ (t ) and p(n,-r) for small t. 

Based on their experimental results, Faureau, Low and 

Pfeffer (1956) hypothesised the distribution of Pq(t) for 

•2 

a gaussian process X(t) whose spectrum is (1 + o 2 ) to 
be negative exponential. However, using his asymptotic ex- 
pression Longuet-Higgins (1962) was able to disprove this 
conjecture. 

Other experimental and analytical approximations of the 
desired probabilities have been given but almost all are 
asymptotic approximations for small r or approximations as 
the level A approaches *. Although, we cannot obtain the 
exact probabilities p(n,x) and p(t), we desire approxima- 
tions which are valid for intermediate level t and n > 1. 
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APPENDIX III 

COMPARISON OF SAMPLE DISTRIBUTIONS 
TO PREDICTED DISTRIBUTIONS 
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■5*5 - •- 


* 


A LEVEL 0.75 AUTOCORRELATION PARAMETER B O.OiOOO 


CUMMULATIVE CUMMULATIVE PREDICTED SAMPLE 


NUMBER OF 

1 

SAMPLE 

I 

PREDICTED 

I 

PROBABILITY 

1 

MINUS 

I 

OVERSHOOTS 

I 

FREQUENCY 

I 

frequency 

I 

U»I) 

I 

PREDICTED 

1 COMMENTS 

0 

■4" 

i 

97.0 

•4" 

I 

73.7 

p 4" 

i 

0.29*8 

• 5 " 

1 

23.27 

I 

1 

i 

136.0 

1 

133.6 

i 

0.2399 

I 

2.30 

I 

2 

i 

163.0 

i 

175.5 

I 

0.1675 

I 

-12.59 

I 

3 

i 

187.0 

i 

203.2 

I 

0.1106 

I 

-16.25 

I 

* 

i 

212.0 

I 

220.9 

I 

0.0709 

I 

-8.98 

1 

5 

i 

226.0 

I 

232.1 

i 

0.0**6 

I 

-6.1* 

I 

fj 

i 

237.0 

I 

239.0 

I 

0.0277 

I 

-2.08 

1 

1 

i 

2*2.0 

I 

2*3.3 

i 

0.0i71 

I 

-1.36 

I 

6 

i 

2*9.0 

I 

2*5.9 

i 

0.010* 

I 

3.02 

1 

9 

i 

250.0 

I 

2*9.9 

l 

0.0169 

I 

0.00 

1 


THE NEGATIVE BINOMIAL MODEL WAS SELECTED. 


ESTIMATED MEAN 
SAMPLE MEAN 


1.9531 

2 . 00*0 


ESTIMATED VARIANCE *.6889 
SAMPLE VARIANCE 5.1365 



A LEVEL 1*50 AUTOCORRELATION PARAMETER B 0*01000 


CUMMULATIVE CUMMUL4T IV.E PREDICTED SAMPLE 


NUMBER OF 
OVERSHOOTS 

I 

I 

SAMPLE 

FREQUENCY 

I 

I 

PREDICTED 

FREQUENCY 

I 

I 

PROBABILITY 

<X*I) 

I 

I 

MINUS 

PREDICTED 

I 

I COMMENTS 

C 

I 

182.0 

I 

182.6 

I 

0.7306 

I 

-0.65 

I 

1 

I 

210.0 

I 

212.6 

I 

0.1199 

I 

-2.64 

I 

2 

I 

230.0 

I 

226.6 

I 

0.0560 

I 

3.33 

I 

3 

I 

236. C 

I 

234.6 

I 

0.0319 

I 

1.36 

I 

4 

I 

240.0 

I 

239.5 

I 

0.0197 

I 

0.42 

I 

5 

I 

242.0 

I 

242.7 

I 

0.0128 

I 

-0.78 

I 

6 

I 

243.9 

I 

244.9 

I 

9.0085 

I 

-1.93 

I 

7 

I 

247.0 

I 

246*4 

I 

0.0058 

1 

0.59 

I 

8 

I 

247.0 

I 

247.4 

I 

0.0040 

I 

-0.42 

I 

9 

I 

250.0 

I 

250.0 

I 

0.0102 

I 

0.00 

I 


THE NEGATIVE BINOMIAL MODEL WAS SELECTED. 


ESTIMATED MEAN 0.7169 
SAMPLE MEAN 0.6920 


ESTIMATED VARIANCE 3.1306 
SAMPLE VARIANCE 2.5352 



A LFVEL 2*00 AUTOCORRELATION PARAMETER B 0*01000 



CUMMULATIVE 

CUMMULATIVE 

PREDICTED 


SAMPLE 

NUMBER OF 

I 

SAMPLE I 

PREDICTED I 

PROBABILITY 

I 

MINUS 

OVERSHOOTS I 

FREQUENCY I 

FREQUENCY I 

<x«n 

I 

PREDICTED 

0 

I 

221.0 I 

216.4 I 

0.8656 

'▼ 

I 

4.57 

1 

I 

231.0 I 

233.6 I 

0.0689 

I 

-2.65 

2 

I 

236.0 I 

240.6 I 

0.0281 

I 

-4.69 

3 

I 

241.0 I 

244.3 I 

0.0145 

I 

-3.33 

4 

I 

243.0 I 

246.4 I 

0.0083 

I 

-3.41 

3 

I 

248.0 I 

247.6 I 

0.0050 

I 

0.31 

6 

I 

249.0 I 

248.4 I 

0.0031 

I 

0.52 

7 

I 

230.0 I 

249.9 I 

0.0061 

I 

0.00 


COMMENTS 


THE NEGATIVE BINOMIAL MODEL WAS SELECTED* 

ESTIMATED MEAN 0.3024 ESTIMATED VARIANCE 1*1484 

SAMPLE MEAN 0.3240 SAMPLE VARIANCE 1.1676 
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( 


A LEVEL 2.50 AUTOCORRELATION PARAMETER B 0.01000 


CUMMULATIVE CUMMULATIVE PREDICTED SAMPLE 


NUMBER OF 
OVERSHOOTS 

I 

I 

SAMPLE 

FREQUENCY 

I 

I 

PREDICTED 

FREQUENCY 

I 

I 

PROBABILITY 
(X»I ) 

I 

I 

MINUS 

PREDICTED 

I 

I COMMENTS 

0 

*♦* 

I 

241.0 

I 

224.1 

I 

0.8966 

I 

16.84 

I 

1 

• 

i 

243.0 

I 

248.6 

I 

0.0978 

I 

-5.62 

I 

2 

I 

246.0 

I 

249.9 

I 

0.0053 

I 

-3.95 

I 

3 

I 

246.0 

I 

250.0 

I 

0.P001 

I 

-4.00 

2 

4 

I 

248.0 

I 

250.0 

I 

0.0000 

I 

-2.00 

I 

5 

I 

248.0 

I 

250.0 

I 

0.0000 

I 

-2.00 

I 

6 

I 

250.0 

I 

250.0 

I 

-0.0000 

I 

0.00 

I 


THE POISSON MODEL WAS SELECTED. 

ESTIMATED MEAN 0.1091 ESTIMATED VARIANCE 0.5158 

SAMPLE MEAN 0.1120 SAMPLE VARIANCE 0.4613 


A LEVEL 3.00 AUTOCORRELATION PARAMETER B 0.01000 


CUMMULATIVE CUMMULATIVE PREDICTED SAMPLE 
NUMBER OF I SAMPLE I PREDICTED I PROBABILITY I MINUS I 

OVERSHOOTS I FREQUENCY I FREQUENCY I (X«I) I PREDICTED I 

0 I 248.0 I 241.7 I 0.9668 I 6.27 I 

1 I 250.0 I 250.0 I 0.0331 I 0.00 I 


THE POISSON MODEL WAS SELECTED. 

ESTIMATED MEAN 0.0336 ESTIMATED VARIANCE 0.2894 

SAMPLE MEAN 0.0080 SAMPLE VARIANCE 0.0079 


COMMENTS 


A LEVEL 0.75 AUTOCORRELATION PARAMETER B 0.05000 


CUMMULATIVE CUMMULATIVE PREDICTED SAMPLE 


NUMBER OF 
OVERSHOOTS 

I 

I 

SAMPLE 

FREQUENCY 

I 

I 

PREDICTED 

frequency 

I 

I 

PROBABILITY 

(Xel) 

I 

I 

MINUS 

PREOICTED 

I 

I COMMENTS 

0 

i 

19.0 

I 

10.8 

“T 

I 

0.0435 

▼ 

I 

8.11 

i 

1 

i 

54.0 

I 

38.7 

I 

0.1115 

I 

15.23 

i 

2 

i 

88.0 

I 

79.2 

I 

0.1619 

I 

8.74 

i 

3 

i 

128.0 

I 

123.0 

I 

0.1751 

I 

4.96 

i 

4 

i 

164.0 

I 

162.2 

I 

0.1570 

I 

1.70 

i 

5 

i 

197.0 

I 

193.1 

I 

0.1233 

I 

3.86 

i 

6 

7 

i 

i 

217.0 

231.0 

• 

I 

215.0 

229.5 

I 

0.0877 

0.0578 

I 

1.92 

1.47 

i 

8 

i 

242.0 

I 

238.* 

I 

0.0357 

I 

3.53 

i 

9 

i 

248.0 

I 

243.7 

I 

0.0210 

I 

4.27 

i 

10 

i 

250.0 

I 

250.0 

I 

0.0250 

I 

0.00 

i 


THE NEGATIVE BINOMIAL MODEL WAS SELECTED. 

ESTIMATED MEAN 3.8901 ESTIMATED VARIANCE 5.9067 

SAMPLE MEAN 3.6480 SAMPLE VARIANCE 5.6266 




A 

LEVEL 1.50 


AUTOCORRELATION PARAMETER 

B 0.05000 




cummulative 

cummulative 

PREDICTED 


SAMPLE 


NUMBER OF 

I 

SAMPLE 

I 

PREDICTED 

I 

PROBABILITY 

I 

MINUS 

I 

OVERSHOOTS 

I 

FREQUENCY 

1 

FREQUENCY 

I 

( X* I ) 

I 

PREDICTED 

I 

0 

▼ 

I 

93.0 

▼ 

I 

84.2 

+ “ 
I 

0.3368 

I 

8.79 

I 

1 

I 

148.0 

I 

144.9 

I 

0.2430 

I 

3.03 

I 

2 

I 

174.0 

I 

184.7 

I 

0.1591 

I 

-10.75 

I 

3 

I 

204.0 

I 

209.9 

I 

0.1006 

I 

-5.92 

I 

4 

I 

230.0 

I 

225.5 

I 

0.0625 

I 

4.43 

I 

5 

I 

235.0 

I 

235.1 

I 

0.0384 

I 

-0.17 

I 

6 

I 

241.0 

I 

241.0 

I 

0.0234 

I 

-0.04 

I 

7 

I 

245.0 

I 

244.6 

I 

0.0142 

I 

0.39 

I 

8 

I 

247.0 

I 

246.7 

I 

0.0086 

I 

0.23 

I 

9 

I 

249.0 

I 

248.0 

I 

0.0051 

I 

0.94 

I 

10 

I 

250.0 

I 

250.0 

I 

0.0077 

I 

0.00 

I 


THE NEGATIVE BINOMIAL MODEL WAS SELECTED* 


ESTIMATED MEAN 1.7514 
SAMPLE MEAN 1.7360 


estimated VARIANCE 
SAMPLE VARIANCE 


4.2506 

4.0906 


COMMENTS 


5 
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A LEVEL 2*00 AUTOCORRELATION PARAMETER B 0.05000 

cummulative cummulative PREDICTED SAMPLE 
NUMBER OF I SAMPLE I °REDICTED I PROBABILITY I MINUS I 

OVERSHOOTS I FREQUENCY I FREQUENCY I (X=I> I PREDICTED I COMMENTS 

— -+ — ♦ — — + — + 


0 

I 

155.0 

I 

140.3 

I 

0.5613 

I 

14.66 

I 

1 

I 

194.0 

I 

199.4 

I 

0.2366 

I 

-5.49 

I 

2 

I 

220.0 

I 

226.3 

I 

0.1074 

I 

-6.35 

I 

3 

I 

236.0 

I 

238.8 

I 

0.0499 

I 

-0.83 

I 

4 

I 

242.0 

I 

244.6 

I 

0.0234 

I 

-2.69 

I 

5 

I 

246.0 

I 

247.4 

I 

0.0111 

I 

-1.47 

I 

6 

I 

249.0 

I 

249.7 

I 

0.0052 

I 

0.20 

l 

7 

I 

250.0 

I 

249.9 

I 

0.0048 

I 

0.00 

I 


THE NEGATIVE BINOMIAL MODEL WAS SELECTED. 

ESTIMATED MEAN 0.8206 ESTIMATED VARIANCE 1.5973 

SAMPLE MEAN 0.8240 SAMPLE VARIANCE 1.8323 


I 


l 

j 
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A LEVEL 2.50 


AUTOCORRELATION PARAMETER B 0.05000 


CUMMULATIVE CUMMULATIVE PREDICTED SAMPLE 


NUMBER OF 
OVERSHOOTS 

I 

I 

sample 

FREQUENCY 

I 

I 

PREDICTED 

FREQUENCY 

I 

I 

PROBABILITY 
(X ;I ) 

I 

I 

MINUS 

PREDICTED 

I 

I 

COMMENTS 

0 

I 

214.0 

m + m 

l 

1 8 1 o 3 

I 

0.7255 

I 

32.61 

I 


1 

I 

238.0 

I 

239.5 

I 

0.2328 

I 

-1.58 

I 


2 

I 

247.0 

I 

248.9 

I 

0.0373 

I 

-1.92 

I 


3 

I 

248.0 

I 

249.9 

I 

0.0039 

I 

-1.92 

I 


4 

I 

248.0 

I 

250.0 

I 

0.0003 

I 

-2.00 

I 


5 

I 

249.0 

I 

250.0 

I 

0.0000 

I 

-1.00 

I 


6 

I 

250.0 

I 

250.0 

I 

-0.0000 

I 

O.CO 

I 



THE POISSON MODEL WAS SELECTED. 




ESTIMATED MEAN 

0.3208 

ESTIMATED 

VARIANCE 

0.7417 

SAMPLE MEAN 

0.2240 

SAMPLE 

VARIANCE 

0.4717 
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A LEVEL 3*00 AUTOCORRELATION PARAMETER B 0.05000 


cummulative cummulative PREDICTED SAMPLE 
NUMBER OF I SAMPLE I PREDICTED I PROBABILITY I MINUS I 

OVERSHOOTS I FREQUENCY I FREQUENCY I (X=I> I PREDICTED I 


0 I 

1 I 

2 I 


241.0 I 

248.0 I 

250.0 I 


225.1 I 
248.7 I 
250.0 I 


0.90C6 I 
0.0942 I 
0.0050 I 


15.84 I 
-0.72 I 
0.00 I 


THE POISSON MODEL WAS SELECTED. 


ESTIMATED MEAN 0.1047 
SAMPLE MEAN 0.0440 


ESTIMATED VARIANCE 0.4330 
SAMPLE VARIANCE 0.0532 


COMMENTS 



A LEVEL 0.75 AUTOCORRELATION PARAMETER 9 0.10000 


CUMMULATIVE CUMMULATIVE PREDICTED SAMPLE 


NUMBER OF 
OVERSHOOTS 

I 

I 

SAMPLE 

FREQUENCY 

I 

I 

PREDICTED 

FREQUENCY 

I 

I 

PROBABILITY 
( X = I ) 

I 

I 

MINUS 

PREDICTED 

I 

I 

COMMENTS 

0 

I 

4.0 

▼ 

I 

1 #6 

▼ 

I 

0.0067 

I 

2.30 

I 


1 

I 

14.0 

I 

9.5 

I 

0.0313 

I 

4.47 

I 


2 

I 

32.0 

I 

28.1 

I 

0.0745 

I 

3.62 

I 


3 

I 

64.0 

I 

58.6 

I 

0.1221 

I 

5.30 

I 


4 

I 

95.0 

I 

97.2 

I 

0.1543 

I 

-2.29 

I 


5 

I 

137.0 

I 

137.4 

I 

0.1607 

I 

-0.47 

I 


6 

I 

170.0 

I 

173.3 

I 

0.1433 

I 

-3.30 

I 


7 

I 

193.0 

I 

201.4 

I 

0.1125 

I 

-8.44 

I 


8 

I 

220.0 

I 

221.3 

I 

0.0794 

I 

-1.30 

I 


9 

I 

237.0 

I 

234.0 

I 

0.0511 

I 

2.91 

I 


10 

I 

245.0 

I 

241.6 

I 

0.0303 

I 

3.32 

I 


11 

I 

248.0 

I 

245.6 

I 

0.0167 

I 

2.12 

I 


12 

I 

249.0 

1 

248.0 

I 

0.0097 

I 

0.94 

I 


13 

I 

250.0 

1 

249.9 

I 

0.0077 

I 

0.00 

I 



THE NEGATIVE BINOMIAL MODEL WAS SELECTED. 


ESTIMATED MEAN 5.4116 
SAMPLE MEAN 5.3680 


ESTIMATED VARIANCE 6.3429 
SAMPLE VARIANCE 6.4584 
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A LEVEL 1*50 AUTOCORRELATION PARAMETER 8 0.10000 


CUMMULATIVE CUMMULATIVE PREDICTED SAMPLE 


NUMBER OP 
OVERSHOOTS 

I 

I 

» 

SAMPLE 

FREQUENCY 

I 

I 

predicted 

FREQUENCY 

I 

I 

PROBABILITY 

(X«I) 

I 

I 

MINUS 

PREDICTED 

I 

I 

COMMENTS 

0 

r 

1 

56.0 

I 

48.9 

T 

I 

0.1956 

1 

9.08 

I 


1 

1 

104.0 

1 

106.6 

I 

0.2310 

I 

-2.68 

I 


2 

I 

142.0 

I 

155.1 

I 

0.1937 

I 

-13.12 

I 


3 

I 

177.0 

I 

190.2 

I 

0. 140 3 

I 

-13.20 

I 


4 

I 

203.0 

I 

213.6 

I 

0.0936 

1 

-10.61 

I 


5 

I 

226.0 

I 

?2e.4 

I 

0.0592 

I 

-2.43 

I 


6 

I 

242.0 

I 

237. u 

I 

0.0361 

T 

4.52 

I 


7 

I 

248.0 

I 

2 42.8 

I 

0.0214 

I 

5.15 

I 


0 

I 

249.0 

I 

245.9 

I 

0.0124 

I 

3.03 

I 


9 

I 

250.0 

I 

250.0 

I 

0.0161 

T 

0.00 

I 



THE NEGATIVE BINOMIAL MODEL WAS SELECTED. 


ESTIMATED MEAN 2.3430 
SAMPLE MEAN 2.4040 


ESTIMATED VARIANCE 4.6488 
SAMPLF VARIANCE 4.3140 



A LEVEL 2.00 AUTOCORRELATION PARAMETER B 0.10000 


CUMMULATtVE CUMMULAT1VE PREDICTED SAMPLE 


NUMBER OP 
OVERSHOOTS 

I 

I 

SAMPLE 

FREQUENCY 

I 

I 

PREDICTED 

FREQUENCY 

I 

I 

PROBABILITY 
(X«I ) 

I 

I 

MINUS 

PREDICTED 

I 

1 

COMMENTS 

0 

I 

124.0 

I 

115.6 

I 

0.4624 

I 

8.38 

I 


1 

I 

179.0 

I 

184.5 

I 

0.2755 

I 

-5.51 

I 


2 

l 

208.0 

I 

219.4 

I 

0.1396 

I 

-11.41 

l 


3 

I 

232.0 

I 

236.0 

I 

0.0665 

I 

-4.05 

I 


A 

I 

243.0 

I 

243.7 

I 

0.0307 

I 

-0.74 

I 


5 

I 

249.0 

I 

247.2 

I 

0.0139 

I 

1.77 

I 


6 

I 

250.0 

I 

250.0 

I 

0.0111 

I 

0.00 

I 



THE NEGATIVE BINOMIAL MODEL WAS SELECTED. 

ESTIMATED MEAN 1.0225 ESTIMATED VARIANCE 1.7547 

SAMPLE MEAN 1.0600 SAMPLE VARIANCE 1.8799 
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A LEVEL 2.50 AUTOCORRELATION PARAMETtR 3 0.10000 

CUMMULATIVE CUMMULATIVE PREDICTED SAMPLE 
NUMBER OF I SAMPLE I PREDICTED I PROBABIL ITY I MINUS I 


OVERSHOOTS 

I 

FREQUENCY 

I 

FREQUENCY 

I 

(X«l ) 

I 

PREDICTED 

I 

0 

I 

195.0 

I 

174.5 

I 

0.6982 

I 

20.44 

* + ■ 
I 

1 

I 

237.0 

I 

237.2 

I 

0.2508 

I 

-0.26 

I 

2 

I 

246.0 

I 

249.5 

2 

0.0450 

I 

-2.52 

I 

3 

I 

249.0 

I 

249.9 

I 

0.0053 

I 

-0.87 

I 

4 

I 

250.0 

I 

249.9 

I 

0.0004 

I 

o.co 

T 

6 


THE POISSON MODEL WAS SELECTED. 

ESTIMATED MEAN 9.3392 ESTIMATED VARIANCE 0.9196 

SAMPLE MEAN 0.2920 SAMPLE VARIANCE 0.4003 


COMMENTS 
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A LEVEL 3*00 AUTOCORRELATION PARAMETER 3 0.10000 


CUMULATIVE CUMMULATIVE PREDICTED SAMPLE 
NUMBER OF I SAMPLE I PREDICTED I PROBABILITY I MINUS I 

overshoots i frequency i frequency i <x-i> i predicted i 


— *■ — - — + — — ♦ — — — — — - — + — + 

0 I 242*0 I 223.8 I 0.9033 I 16.15 I 

1 I 249.0 I 2<*8.7 I 0.0917 I 0.20 I 

2 » 250.0 I 249.9 1 0.0048 I 0.00 I 


THE POISSON MODEL WAS SELECTED. 

ESTIMATED MEAN 0.1016 ESTIMATED VARIANCE 0.**813 

SAMPLE MEAN 0.0360 SAMPLE VARIANCE 0.0428 
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COMMENTS 



A LEVEL 0*75 AUTOCORRELATION PARAMETER B 0.5000C 


CUMMULATlvE CUMmulATIvE PREDICTED SAMPLE 


number of 

OVERSHOOTS 

I 

I 

SAMPLE 

FREQUENCY 

I 

I 

PREDICTED 

FREQUENCY 

I 

I 

PROBABILITY 

(X*I) 

I 

I 

MINUS 

PREDICTED 

I 

I 

COMMENTS 

0 

I 

0.0 

i 

0.0 

m + m 
I 

0.0000 

I 

-0.00 

I 


1 

I 

0.0 

i 

0.0 

I 

0.0000 

I 

-0.01 

I 


2 

1 

0.0 

i 

0.0 

I 

0.0003 

I. 

-0.09 

I 


3 

I 

0.0 

i 

0.4 

I 

0.0015 

I 

-0.48 

I 


4 

I 

0.0 

i 

1.8 

I 

0.0056 

I 

-1.89 

I 


5 

I 

2.0 

i 

5.7 

I 

0.0155 

I 

-3.77 

I 


6 

I 

9.0 

i 

14.3 

I 

0.0343 

I 

-5.35 

T 

4 


7 

I 

26.0 

i 

30.0 

I 

0.0628 

I 

-4.07 

I 


8 

I 

53.0 

l 

54.3 

I 

0.096 9 

I 

-1.32 

I 


9 

I 

77.0 

i 

86.2 

I 

0.1278 

I 

-4.29 

I 


10 

I 

111.0 

i 

122.7 

I 

0.1456 

I 

-li.71 

I 


11 

I 

1*9.0 

i 

15*-. 8 

I 

0.1445 

I 

-9.86 

I 


12 

I 

173.0 

i 

190.3 

I 

0.1255 

I 

-17.31 

T 

A 


13 

I 

202.0 

i 

214.4 

I 

0.0964 

I 

-12.43 

I 


u 

I 

222.0 

i 

230.7 

I 

0.0654 

I 

-8.79 

I 


15 

I 

235.0 

I 

240.6 

I 

0.0393 

I 

-5.62 

I 


16 

I 

244.0 

i 

245.8 

I 

0.021C 

I 

-1.37 

I 


17 

I 

248.0 

i 

248.3 

I 

0.0099 

I 

-0.36 

* 

1 


18 

I 

249.0 

l 

2*9.4 

I 

0.0042 

I 

-0.42 

I 


19 

I 

249.0 

l 

249.8 

I 

0.C015 

I 

-0.81 

I 


20 

I 

249.0 

i 

249.9 

I 

0.0005 

I 

-0.94 

I 


21 

I 

250.0 

i 

250.0 

I 

0.0002 

T 

C.00 

I 



THE BINOMIAL MODEL WAS SELECTED. 

ESTIMATED MEAN 10. ,697 ESTIMATED VARIANCE 7.2353 

SAMPLE MEAN 11.0030 SAMPLE VARIANCE 7.7830 



J J 



A LEVEL 1*50 AUTOCORRELATION PARAMETER B 0.50000 


CUMMULATIVE CUMMULATIVE PREDICTED SAMPLE 


NUMBER OF 
OVERSHOOTS 

I SAMPLE 
I FREQUENCY 

I 

I 

PREDICTED 

FREQUENCY 

I 

I 

PROBABILITY 
(X*I J 

I 

I 

MINUS 

PREDICTED 

I 

I 

COMMENTS 

0 

I 3.0 

T 

I 

5.8 

I 

0.0233 

i 

-2.83 

I 


1 

I 7.0 

I 

25.2 

I 

0.0777 

i 

-18.26 

I 


2 

I 35.0 

I 

59.7 

I 

0.1378 

i 

-24.72 

I 


3 

I 84.0 

I 

102.9 

I 

0.1730 

i 

-18.99 

I 


4 

I 129.0 

l 

146.1 

I 

0.1724 

i 

-17.12 

I 


5 

I 172.0 

I 

182.3 

I 

0.1450 

i 

-10.38 

I 


6 

I 205.0 

I 

209.1 

I 

0.1069 

i 

-4.13 

I 


7 

I 229.0 

I 

226.8 

I 

0.0709 

i 

2.12 

I 


8 

I 241.0 

I 

237.6 

I 

0.0431 

i 

3.33 

I 


9 

I 247.0 

I 

243.7 

I 

0.0243 

I 

3.24 

I 


10 

I 247.0 

I 

246.9 

I 

0.0129 

i 

0.01 

I 


11 

I 249.0 

I 

248.6 

I 

0.0064 

i 

0.39 

I 


12 

I 250.0 

I 

250.0 

I 

0.0055 

i 

0.00 

I 



THE NEGATIVE BINOMIAL MODEL WAS SELECTED. 

ESTIMATED MEAN 4.2630 ESTIMATED VARIANCE 5.4600 

SAMPLE MEAN 4.6080 SAMPLE VARIANCE 4.2071 


\ •• 


A LEVEL 2*00 AUTOCORRELATION PARAMETER B 0.50000 


CUMMULATIVE CUMMULATIvE PREDICTED SAMPLE 


NUM8ER OF 
OVERSHOOTS 

I 

I 

SAMPLE 

FREQUENCY 

I 

I 

PREDICTED 

FREQUENCY 

I 

I 

PROBABIL ITY 
< X=I > 

I 

I 

MINUS 

PREDICTED 

I 

I 

COMMENTS 

0 

I 

50.0 

*+' 

I 

59.0 

I 

0.2362 

I 

-9.06 

I 


1 

I 

128*0 

I 

134.8 

I 

0.3030 

I 

-6.82 

I 


2 

I 

183.0 

I 

191.5 

I 

0.2267 

I 

-8.50 

I 


3 

I 

226.0 

I 

223.6 

I 

0.1291 

I 

2.19 

I 


6 

I 

242.0 

I 

239.3 

I 

0.0621 

I 

2.67 

I 


5 

I 

247.0 

I 

2 45.9 

I 

0.0265 

I 

1.03 

I 


6 

I 

248.0 

I 

248.5 

I 

0.0103 

I 

-0.56 

I 


7 

I 

249.0 

I 

249.5 

I 

0.0038 

I 

-0.51 

I 


8 

I 

250.0 

I 

250.0 

I 

0.0019 

I 

-0.00 

I 



THE NEGATIVE BINOMIAL 


MODEL WAS 


SELECTED. 


ESTIMATED MEAN 1.6308 
SAMPLE MEAN 1.7080 


ESTIMATED VARIANCE 2.0733 
SAMPLE VARIANCE 1.9505 



A LEVEL 2*50 AUTOCORRELATION PARAMETER B 0.50000 


CUMMULATIVE CUMMULATtVE PREDICTED SAMPLE 


NUMBER OF 

I 

SAMPLE 

I 

PREDICTED 

I 

PROBABILITY 

I 

MINUS 

OVERSHOOTS 

I 

A t 

FREQUENCY 

I 

FREQUENCY 

I 

U«I) 

I 

PREDICT 

0 

T 

I 

155.0 

I 155.9 

I 

0.6239 

I 

-0.97 

1 

I 

226.0 

I 

229.5 

I 

0.2943 

I 

-3.56 

2 

I 

246.0 

I 

246.9 

I 

0.0694 

I 

1.07 

3 

I 

250.0 

I 

249.9 

1 

0.0123 

I 

0.00 


the poisson model was selected. 

ESTIMATED MEAN 0.4717 ESTIMATED VARIANCE 0.9764 

SAMPLE MEAN 0.4640 SAMPLE VARIANCE 0.4756 


COMMENTS 
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A LEVEL 3*00 AUTOCORRELATION PARAMETER 9 0.50000 


CUMULATIVE CUMMULATIVE PREDICTED SAMPLE 
NUMBER OF I SAMPLE I PREDICTED I PROBABILITY I MINUS I 

OVERSHOOTS I FREQUENCY I FREQUENCY I (X«I) I PREDICTED I 


---- + 

0 I 

1 I 

2 I 


222*0 I 
248*0 I 
250*0 I 


225*4 I 
248*7 I 
249*9 I 


0*9019 I 
0*0930 I 
0*0049 I 


-3*49 I 
-0.76 I 
0.00 I 


THE POISSON MODEL WAS SELECTED* 


ESTIMATED MEAN 
SAMPLE MEAN 


0*1032 

0.1200 


ESTIMATED VARIANCE 0.5785 
SAMPLE VARIANCE 0.1220 


COMMENTS 
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A LEVEL 0*75 AUTOCORRELATION PARAMETER B 1.00000 


CUMMULATIVE CUMMULATIVE PREDICTED SAMPLE 


NUMBER OF 
OVERSHOOTS 

I 

I 

SAMPLE 

FREQUENCY 

I 

I 

PREDICTED 

FREQUENCY 

I 

I 

PROBABILITY 

(X*I) 

I 

I 

MINUS 

PREDICTED 

I 

I COMMENTS 

0 

"♦* 

T 

A 

0.0 

I 

0.0 

I 

0.0000 


-0.00 

' 

► 

1 

I 

0.0 

I 

0.0 

I 

0.0000 

i 

-0.00 

I 

2 

I 

0.0 

I 

0.0 

I 

0.0000 

i 

-0.00 

I 

3 

I 

0.0 

I 

0.0 

I 

0.0000 

i 

-0.00 

I 

4 

I 

0.0 

I 

0.0 

1 

0.0000 

i 

-0.02 

I 

5 

I 

0.0 

I 

0.1 

I 

0.0004 

i 

-0.13 

I 

6 

I 

1.0 

I 

0.5 

I 

0.0015 

i 

0.46 

I 

7 

I 

2.0 

I 

1.7 

I 

0.0047 

i 

0.26 

I 

• 

I 

5.0 

I 

4.7 

I 

0.0120 

i 

0.24 

I 

9 

I 

13.0 

I 

11.2 

I 

0.0259 

i 

1.76 

I 

10 

I 

24.0 

I 

23.2 

I 

0.0478 

i 

0.79 

I 

11 

I 

41.0 

I 

42.4 

I 

0.0768 

i 

-1.40 

I 

12 

I 

72.0 

I 

69.2 

I 

0.1075 

i 

2.70 

I 

13 

I 

100.0 

1 

102.2 

I 

0.1320 

i 

-2.29 

I 

14 

I 

132.0 

I 

137.9 

I 

0.1425 

i 

-5.94 

I 

15 

I 

172.0 

I 

171.8 

I 

0.1357 

i 

0.11 

I 

16 

I 

195.0 

I 

200.3 

I 

0.1140 

i 

-5.39 

I 

17 

I 

225.0 

I 

221.5 

I 

0.0845 

i 

3.47 

I 

18 

I 

234.0 

I 

235.3 

1 

0.0552 

i 

-1.32 

I 

19 

I 

239.0 

I 

243.2 

I 

0.0317 

i 

-4.26 

I 

20 

I 

242.0 

I 

247.2 

I 

0.0159 

i 

-3.26 

I 

21 

I 

249.0 

I 

249.0 

I 

0.0070 

i 

-0.02 

I 

22 

I 

250.0 

I 

250.0 

i 

0.0039 

i 

-0.00 

I 


THE BINOMIAL MODEL WAS SELECTED. 

ESTIMATED MEAN 13.9498 ESTIMATED VARIANCE 7.5813 

SAMPLE MEAN 14.2160 SAMPLE VARIANCE 4.3386 
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A LEVEL 1.50 AUTOCORRELATION PARAMETER B 1.00000 


CUMMULATIVE CUMMULATIVE PREDICTED SAMPLE 


NUMBER OF 
OVERSHOOTS 

I 

I 

SAMPLE 

FREQUENCY 

I 

I 

PREDICTED 

FREQUENCY 

I 

I 

PROBABILITY 
( X = I ) 

I 

I 

MINUS 

PREDICTED 

I 

I 

COMMENTS 

0 

•4- * 

i 

0.0 

I 

1.3 

“4* 

i 

0.0053 

“T* 

I 

-1.33 

I 


1 

i 

6#0 

I 

8.0 

i 

0.0270 

I 

-2.08 

I 


2 

I 

21.0 

I 

25.4 

i 

0.0693 

I 

-4.41 

I 


3 

i 

46.0 

I 

55.4 

i 

0.1199 

I 

-9.40 

I 


4 

i 

73.0 

I 

94.3 

i 

0.1576 

I 

-21.81 

I 


5 

i 

118.0 

I 

136.7 

i 

0.1677 

I 

-18.74 

I 


6 

i 

169.0 

I 

174.3 

i 

0.1504 

I 

-5.35 

I 


7 

i 

199.0 

I 

203.6 

i 

0.1170 

I 

-4.61 

I 


8 

i 

22<*.0 

I 

223.7 

i 

0.0606 

I 

0.23 

I 


9 

i 

238.0 

I 

236.2 

i 

0.0499 

I 

1.75 

I 


10 

i 

244.0 

I 

243.2 

i 

0.0281 

I 

0.72 

I 


11 

i 

248.0 

I 

246.9 

i 

0.0145 

I 

1.08 

I 


12 

i 

249.0 

I 

248.6 

i 

0.0069 

I 

0.33 

I 


13 

i 

250.0 

I 

249.9 

i 

0.0053 

I 

0.00 

I 



THE NEGATIVE BINOMIAL MODEL WAS SELECTED. 


ESTIMATED MEAN 5. *034 
SAMPLE MEAN 5.6600 


ESTIMATED VARIANCE 
SAMPLE VARIANCF 


5.7737 

5.1248 



A LEVEL 2*00 AUTOCORRELATION PARAMETER B 1*00000 


CUMMULATIVE CUMMULATIVE PREDICTED SAMPLE 


NUMBER OF 
OVERSHOOTS 

I 

I 

SAMPLE 

FREQUENCY 

I 

I 

PREDICTED 

FREQUENCY 

1 

I 

PROBABILITY 

(X*I) 

I 

I 

MINUS 

PREDICTED 

I 

I COMMENTS 

0 

I 

33*0 

I 

38*5 

m + m 
I 

0*1543 

T 

I 

-5.59 

I 

1 

I 

103*0 

I 

106*9 

I 

0*2734 

I 

-3*96 

I 

2 

I 

156*0 

I 

171*0 

I 

0*2561 

I 

-15.00 

I 

3 

I 

206*0 

I 

213*1 

I 

0*1686 

I 

-5.17 

I 

4 

I 

238*0 

I 

235*0 

I 

0*0875 

I 

2*93 

I 

5 

I 

244*0 

I 

244*6 

I 

0*0381 

I 

-0.6C 

I 

6 

I 

249*0 

I 

248*2 

I 

0*0145 

I 

0*76 

I 

7 

I 

250*0 

I 

249*9 

I 

0*0070 

I 

0*00 

I 


THE NEGATIVE BINOMIAL MODEL WAS SELECTED. 


ESTIMATED MEAN 1*9724 
SAMPLE MEAN 2*0760 


ESTIMATED VARIANCE 2*1960 
SAMPLE VARIANCE 2*1347 


A LEVEL 0*75 AUTOCORRELATION PARAMETER B 3*00000 


CUMMULATIVE CUMULATIVE 


NUMBER OF 
OVERSHOOTS 

I 

I 

SAMPLE 

FREQUENCY 

I 

I 

PREDICTED 

FREQUENCY 


0 

I 

0.0 

m + m 
I 

0.0 


1 

I 

0.0 

I 

0.0 


2 

I 

0.0 

I 

0.0 


3 

l 

0.0 

I 

0.0 


4 

I 

0.0 

I 

0.0 


5 

I 

0.0 

I 

0.0 


6 

I 

0.0 

I 

0.0 


7 

I 

0.0 

I 

0.0 


8 

I 

0.0 

I 

0,1 


9 

I 

0.0 

I 

0.4 


10 

I 

1.0 

I 

1.3 


11 

I 

2.0 

I 

3.5 


12 

I 

12.0 

I 

e.i 


13 

I 

23*0 

I 

16.9 


14 

I 

36.0 

I 

31.6 


15 

I 

59.0 

I 

53.2 


16 

I 

82.0 

I 

81.7 


17 

I 

111.0 

I 

114.9 


18 

I 

159.0 

I 

149.4 


19 

I 

193*0 

I 

181.2 


20 

I 

207.0 

I 

20Y.2 


21 

I 

221.0 

I 

226.0 


22 

I 

237.0 

I 

238.0 


23 

I 

244.0 

I 

244.6 


24 

I 

246.0 

I 

247.9 


25 

I 

249.0 

I 

249.3 


26 

I 

250.0 

I 

250.0 



THE BINOMIAL MOOEL WAS SELECTED. 

ESTIMATED MEAN 17.5446 
SAMPLE MEAN 17*6720 


PREDICTED SAMPLE 
PROBABILITY I MINUS I 

( X* I ) I PREDICTED I COMMENTS 


0*0000 I -0*00 I 
0*0000 I -0*00 I 
0*0000 I -0*00 I 
0*0000 I -0*00 I 
0.0000 I -0.00 I 
0.0000 I -0.00 I 
0.0000 I -0.00 I 
0.0001 1 -0.03 1 
0.0003 I -0.13 I 
0.0012 I -0.44 I 
0.0035 I -0.34 I 
0.0087 I -1.51 I 
0.0186 I 3.82 I 
0.0351 1 6.03 I 
0.0586 I 4.39 I 
0.0866 I 5.71 I 
0.1138 I 0.24 I 
0.1329 I -3.99 I 
0.1379 I 9.51 I 
0.1271 I 11.71 I 
0.1039 I -0.27 I 
0.0751 I -5.05 I 
0.0478 I -1.01 I 
0.0267 I -0.69 I 
0.0129 I -1.94 I 
0.0054 I -0.31 I 
0.0027 I 0.00 I 


VARIANCE 8.0941 
VARIANCE 8.7935 


estimated 

SAMPLE 


A LEVEL 1*50 AUTOCORRELATION PARAMETER B 3.00000 


CUMMULATIVE CUMMULATIVE PREDICTED SAMPLE 


NUMBER OF 
OVERSHOOTS 

I 

I 

SAMPLE 

FREQUENCY 

I 

I 

PREDICTED 

FREQUENCY 

1 

I 

PROBABILITY 

<X*I) 

I 

I 

MINUS 

PREDICTED 

I 

1 

COMMENTS 

0 

▼ 

I 

0.0 

I 

0.3 

► — 

0.0015 

• + * 
I 

-0.37 

I 


1 

I 

1.0 

I 

2.8 

I 

0.0098 

I 

-1.84 

I 


2 

I 

6#0 

I 

10.8 

I 

0.0319 

I 

-4.83 

I 


3 

2 

24.0 

I 

28.1 

I 

0.0691 

I 

-4.13 

♦ 

l 


4 

I 

49.0 

I 

56.2 

i 

0.1122 

I 

-7.20 

I 


5 

I 

87.0 

I 

92.6 

I 

0.1457 

I 

-5.62 

I 


6 

I 

140.0 

I 

132.0 

I 

0.1576 

I 

7.96 

I 


7 

I 

175.0 

I 

168.5 

I 

0.1461 

I 

6.43 

I 


8 

I 

209.0 

I 

198.1 

I 

0. I 1 P 5 

I 

10.80 

I 


9 

I 

232.0 

I 

219.5 

I 

0.0854 

I 

12.43 

I 


10 

I 

241.0 

I 

233.4 

I 

0.0554 

I 

7.56 

I 


11 

I 

247.0 

I 

241.6 

I 

0.032 7 

I 

5.38 

I 


12 

I 

248.0 

I 

246.0 

I 

0.0177 

I 

1.95 

I 


13 

I 

249.0 

I 

248.2 

I 

0.0088 

I 

0.74 

I 


14 

I 

250.0 

I 

250.0 

I 

0.0069 

I 

o.oc 

I 



THE POISSON MODEL WAS SELECTED. 

ESTIMATED MEAN 6.4899 ESTIMATED VARIANCE 6.2380 

SAMPLE MEAN 6.3680 SAMPLE VARIANCE <*.7877 
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I 


A LEVEL 2*00 AUTOCORRELATION PARAMETER B 3.00000 


CUMMULATIVE CUMMuLATIvE predicted sample 


NUMBER OF 
OVERSHOOTS 

I 

I 

SAMPLE 

FREQUENCY 

I 

I 

PREDICTED 

FREQUENCY 

I 

I 

PROBABILITY 

(X»!) 

I 

I 

MINUS 

PREDICTED 

I 

I 

COMMENTS 

0 

I 

23*0 

I 

22*6 

I 

0*0906 

I 

0*34 

T 

I 


1 

I 

79*0 

I 

77.0 

I 

0*2175 

I 

1*95 

I 


2 

I 

137.0 

I 

142*3 

I 

0*2612 

I 

-5.34 

I 


3 

I 

193*0 

I 

194*6 

I 

0*2090 

I 

-1*61 

I 


4 

I 

230*0 

I 

225*9 

I 

0*1255 

I 

4. CO 

I 


5 

I 

242*0 

I 

241.0 

I 

0*0602 

I 

0*93 

I 


6 

I 

249*0 

I 

247.0 

I 

0*0241 

I 

1*90 

I 


7 

I 

250*0 

I 

249.9 

I 

0*0116 

I 

0*00 

I 



THE POISSON MODEL WAS SELECTED* 

ESTIMATED MEAN 2*4012 ESTIMATED VARIANCE 2.3773 

SAMPLE MEAN 2*3860 SAMPLE VARIANCE 2*2464 




A LEVEL 2.50 AUTOCORRELATION PARAMETER B 3.00000 





CUMMULATIVE 


CUMMULATIVE 

PREDICTED 


SAMPLE 


NUMBER OE 

I 

SAMPLE 

I 

RREDICTED I 

PROBABILITY 

I 

MINUS 

I 

OVERSHOOTS 

I 

FREQUENCY 

I 

+ 

FREQUENCY I 

(X»I ) 

I 

PREDICTED 

I 

0 

I 

146*0 

I 

126*4 I 

0.5057 

I 

17.55 

I 

1 

I 

226*0 

I 

212*6 I 

0.3447 

I 

13*35 

I 

2 

I 

246*0 

I 

242*0 I 

0.1175 

I 

3*98 

I 

3 

I 

250*0 

I 

*49. 9 I 

0*0319 

I 

0*00 

I 


COMMENTS 


THE POISSON MODEL WAS SELECTED* 

ESTIMATED MEAN 0*6816 ESTIMATED VARIANCE 1*1250 

SAMPLE MEAN 0*5360 SAMPLE VARIANCE 0.5067 




A LEVEL 3.00 AUTOCORRELATION PARAN _ TER B 3.00000 


CUMULATIVE CUMMULATIVE PREDICTED SAMPLE 
NUMBER OF I SAMPLE I PREDICTED I PROBABILITY I MINUS I 

OVERSHOOTS I FREQUENCY I FREQUENCY I ( X- I » I PREDICTED I 

— — ——————— — * *♦- ♦ 

0 1 236*0 I 215.5 I 0.8620 I 20.48 I 

1 I 250.0 I 250.0 I 0.1379 I 0.00 I 


THE POISSON MODEL WAS SELECTED* 

ESTIMATED MEAN 0.1484 ESTIMATED VARIANCE 0.6699 

SAMPLE MEAN 0*0560 SAMPLE VARIANCE 0*0530 


COMMENTS 
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